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ON THE POPULATIONS OF RADIO GALAXIES WITH EXTENDED MORPHOLOGY AT z < 0.3 
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ABSTRACT 

Extended extragalactic radio sources have traditionally been classified into 'Fanaroff & RilevI (FR) I and II 
types, based on the ratio r, of the separation S between the brightest regions on either sides of the host galaxy 
and the total size T of the radio source (r, = S/T). In this paper we examine the distribution of various 
physical properties as a function of of 1,040 luminous (L > L^,) extended radio galaxies (RGs) at z < 0.3 
selected with well-defined criteria from the SDSS, NVSS, and FIRST surveys. About 2/3 of the RGs are 
lobe-dominated (LD), and 1/3 have prominent jets. If we follow the original definition of the FR types, i.e., 
a division based solely on r,, FR I and FR II RGs overlap in their host galaxy properties. However, the 
rare, LD sources with > 0.8 and [O III]A5007 line luminosity > lO^L© are markedly different on average 
from the rest of the RGs, in the sense that they are hosted in lower mass galaxies, live in relatively sparse 
environments, and likely have higher accretion rates onto the central supermassive black hole (SMBH). Thus 
these high emission line luminosity, high LD RGs and the rest of RGs form a well-defined dichotomy. 
Motivated by the stark differences in the nuclear emission line properties of the RG subsamples, we suggest 
that the accretion rate onto the SMBH may play the primary role in creating the different morphologies. At 
relatively high accretion rates, the accretion system may produce powerful jets that create the "classical double" 
morphology (roughly corresponding to the LD sources with > 0.8 and emission lines); at lower accretion 
rates the jets from a radiatively inefficient accretion flow generate radio lobes without apparent "hot spots" at 
the edge (corresponding to the majority of LD sources). At slightly lower accretion rates and in galaxies with 
dense galactic structure, sources with prominent jets result. It is possible that while the high accretion rate 
systems could affect sub-Mpc scale environments, the jets from lower accretion rate systems may efficiently 
suppress activity within the host galaxies. 

Subject headings: galaxies: active - galaxies: elliptical and lenticular, cD - radio continuum: galaxies 
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1. INTRODUCTION 



Ever since the seminal work by iFanaroff & RilevI (11974 
hereafter FR), radio galaxies (RGs) have been broadly catego- 
rized into two types according to their large-scale radio mor- 
phology (e.g., De Y oungll2602) . FR proposed a simple binary 
classification scheme, based on the ratio rj of the separation 
between the brightest regions on either sides of the RG and 
the total size of the radio source. If rj > 0.5, the source was 
considered type II ("edge-brightened"); otherwise the source 
belonged to type I ("ed ge-darkened"). FR found that nearly 
all 3CR (Macka3ll9lil) radio sources with power (P) at 178 
MHz greater than 1 .3 x 10^^ W/Hz were type II, while nearly 
all of those weaker than this power were type I. 

Subsequent studies have found an overlap in radio power 
(about two orders of magnitude) for the transition from one 
type to the other (e.g., Baum & Heckman 1989). In a se- 
ries of papers, Owen and collaborators suggested this divi- 
sion between type I and II RGs was a function of the optical 
luminosity (L) of the host galaxies, in the sense that the di- 
vision was at higher radio power for more luminous galaxies 
{P oc L^MO wen & Laing 1989; Owen & White 1991; Owe^ 
119931; iLedlow & Owenl ll996). A coroflai-y is that, at fixed 
radio power, type I sources are hosted in more optically lumi- 
nous galaxies than the type lis. At a given optical luminosity, 
Owen et al. suggested that the transition as a function of radio 
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power was quite abrupt. 

It was also found that the two types of sources exhibit 
differences other than the radio power. Type lis are usu- 
ally found in less dense environments, and strong emis- 
sion lines c an often be seen in their optic al spectra (e 
IZirbel & Ba um 1995; Zirbel 1997; KauffmannetiD llQOi 
type Is, on the other hand, are usually hosted by giant el- 
liptical galaxies, and on average have weaker or no opti- 
cal nuclear emission lines, which puts them on a differ- 
ent radio pow er-emission line luminosity correlation than 
that of FR lis jHine & Longaid [19791; IZirbel & BaumI [19951 
see below). The cosmological evolution of the two types 
may also differ sig nificantly (e.g., Willott et al. 2001; but see 
iGendre et alJl2010h . Collectively, these differences are usu- 
ally referred to as the FR I/II dichotomy, which has been 
the focus of num erous studies (e.g., [H eckman et al. 1994 
'Baum et al.' '1995; 'Hardcastle etalj[2007i; iKauffmann et atl 
2008; Baldi & Capetti 2010, and references therein). 

Larger samples and better observations have led to sev- 
eral proposals for refinement/alternatives of radio sourc e clas- 
sification sc hemes (e.g., O wen &^_Lsingj 1989; Leahy 119931 : 
lLainglll993l) . In particular, Owen & LaingI (11989 . hereafter 
OL89) categorized RGs into three types based on their radio 
morphologies: classical double (CD) sources are those with 
"compact outer hotspots and elongated, diffuse lobes extend- 
ing from the hotspots back toward the nucleus"; fat double 
(FD) objects have "bright outer rims of radio emission and 
roundish diffuse radio lobes"; while twin jet (TJ) RGs "can 
be described by symmetric jets originating in the nucleus and 
extending on both sides" of the host. OL89 suggested that the 
hosts of FD and TJ sources have similar optical properties. 
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and regarded both types, together with the Narrow-Ang le Tail 
(NAT) objects, as FR I (see also lOwen & Whitelll99ll . The 
CD objects were considered to be equivalent to FR II. 

In addition to the radio morphologies, RGs have been clas- 
sified based on their nuclear optical narrow em ission line 
properties (e.g., i Laing et all Il994t |Ha rdcastle" et al. 200^. 
Objects with weak emission lines were generally referred to 
as low-excitation (LE) RGs, while their counterparts with 
strong lines were known as high-excitation (HE) RGs. Such 
a scheme is believed to better reflect any differences in the 
central supermassive black hole (SMBH), or/and the physi- 
cal c onditions of the a ccretion flow onto the central engine 
(e.g.. iKauffmann et al.l 12008) . On the other hand, classifi- 
cation based on radio morphology is Ukely more intimately 
connected to interactions of the jets with the ( large-scale) 
environments (e.g., De Young 1993; K awakatu e t al. 2009). 
As such, the correspondence between HE/LE and FR II/I 
is not perf ect; a large fractio n of FR lis have a LE nu- 
cleus (e.g .. ILaing et al.l 119941). while some FR Is are HE 
RGs (see iHevwood et al.l l2007h . This suggests that a hy- 
brid classification system that incorporates both radio mor- 
phology and nuclear emission line activity may perform bet- 
ter in revealing distinct populations of RGs, which in turn 
would lead to a fuller understanding of the onset of radio ac- 
tivity, as well as the unificati on of radio-loud a ctive galac- 
tic nuclei (AGNs; e g. , BartheJl98^JUrrv & Pad ovaniiri99"l 
iFalcke & Biermannlll995l; iHardcastle et al.i i2006). 

The samples used by many of the previous studies on 
the FR I/II dichotomy were limited by the available data at 
the time, and thus often have heterogeneous or igins and did 
not have well-define d selection criteria (e.g., lOwenI 119931; 
IZrrbel & BaumlfT995t) . Furthermore, it is inevitable that clas- 
sification of extended radio sources will be subject to the clas- 
sifiers' experience and preference. Both factors may make it 
difficult to uncover the origin of the dichotomy. 

Our primary goal is to understand the origin of the differ- 
ent radio morphologies, which drives us to investigate ways 
to separate RGs into distinct populations, by exploring vari- 
ous classification schemes based on radio morphology or/and 
nuclear activity. In addition, we will use rj as a continuous 
parametrization of the radio morphology to study the transi- 
tion from FR I to FR II, and to quantify any bimodality in the 
physical properties of the host galaxies on the radio power- 
optical luminosity plane. We will mainly use our own termi- 
nology to refer to the RG populations identified in this paper; 
unless specifically noted, we adopt the original definition of 
FR when we refer to the FR I/II types of extended RGs. 

We have constructed a large RG sample with well-defined 
criteria (see and paid particular attention to the measure- 
ments of total size {T) of a RG, as well as the separation (S) 
between the highest surface brightness spots (©. From these 
quantities we define the ratio objectively and reproducibly, 
and present a classification system that is purely based on ra- 
dio morphologies (©. In ^we examine the distribution of 
RGs on the radio power-optical luminosity plane as a function 
of Tj, while in ^we study dependences of various physical 
properties of host galaxies as a function of r,. Comparisons 
between radio loud (RL) and radio quiet (RQi3 galaxies are 
made in ^5.41 Discussion of the FR dichotomy, a proposal for 



an improved classification scheme based on both and nu- 
clear emission line strengths, comparison of our scheme with 
those of FR and OL89, and the possible origin of the mor- 
phological differences, are presented in ^ We summarize 
our main results in ^ Possible systematics due to the limi- 
tations of the radio data we use, and the uncertainties in the 
measurement of r,, are discussed in the Appendices. 

Most of the radio sources in the local Universe are com- 
pact or barely resolved at ~ 5" resolution (see ©; this paper 
concentrates on extended radio sources, and we will study the 
physical properties of the host galaxies of compact sources in 
a future publication. 

Throughout this paper we adopt a flat ACDM cosmological 
model where VIm = 1 - i^A = 0.3 and Ho = 100/ikms~'Mpc~' 
with h = 0.70. As our sources are at z < 0.3, the small dif- 
ferences between the ado pted cosmology and th e current best 
concordance model (e.g.. !Komatsu et ani2010l) have a negli- 
gible effect on our results. 

2. THE PARENT RADIO GALAXY SAMPLE 

We construct the RG sam ple by cros s-matching the Sloan 
Digital Sky Survey (S DSS; lYorket al.1 2000) main galaxy 
spectroscopic sample (!Strauss et al.l 2002) with the NRAO 
VLA Sky Survey (NVSS; Condo n et al.ni998l) and FIRST 

tint Images of the Radio Sky at Twenty-Centimeters; 
:keretal. 1995) surveys. We start with the DR6 
lelman-McCarthy etal.il2008h version of flie NYU Value- 
Added Galaxy CatalosH (VAGC; iBlanton et al.ll200 5l). select- 
ing from the large-scale structure bbrightO subsampl^U 
229,379 galaxies with 0.02 < z < 0.3 and ' < -21 .27 (i.e., 
more luminous than the ch aracteristic magni tude M* in the 
galaxy luminosity function, Blanton et al. 2003. Mj?-' or r" ' 
denotes the SDSS r-band shifted blueward by a factor of 1.1 
in wavelength. Similarly, ' and ' refer to the shifted g and 
/ bands, respectively). Selecting RGs with a uniform absolute 
magnitude limit makes it easy to compute the fraction of RL 
objects in volume-limited galaxy samples. Since the parent 
VAGC galaxy sample excludes quasars, our sample does not 
contain radio quasars. Our survey area covers 6008 deg^. 

As the first step, automatic matching to the ra dio source 
catalog s is carried out following the prescription of lBest et alJ 
( l2005al) . In short, if a galaxy has only one NVSS source pro- 
jected within 3', the pair would be matched depending on 
their angular separation as well as the properties of the FIRST 
source(s) (if present) in the vicinity of the galaxy. If there are 
at least two NVSS sources within 3' of a galaxy, the match- 
ing depends on the spatial distribution and the fluxes of the 
NVSS and FIRST sources. We then visually inspect all galax- 
ies with at least one NVSS source within 3' (irrespective of 
the results of auto-matching), correcting for any mis-match, 
keeping only RGs with total flux density /i,4 > 3 mJy at 1.4 
GHz, and recording their morphological and structural infor- 
mation wherever possible (see The radio flux from the 
NVSS catalog is in general adopted, as fluxes from extended 
sources may be resolved out by the FIRST survey. For com- 
plicated sources that are blended/unresolved in the NVSS im- 
ages, we make use of both NVSS and FIRST data to assign 
proper fluxes to individual RGs. More details on the construc- 
tion of the RG sample will be presented in a future publication 



' We do not use an optical-to-radio luminosity ratio to distinguish between 
RL and RQ galaxies; by RL we refer to galaxies with 1.4 GHz power P > 
10^3 W/Hz, while byRQ we mean galaxies not detected by the NVSS survey 
(i.e., 20cm flux density < 3 mJy; see 0. 



^'http : / / sdss . physics . nyu . edu/ vagc/ 

^ T his subsample is defined by a constant, extinction-coiTected IPetrosiaJ 



(T923) magnitude limit of r = 17.6, together with no fiber collision correc- 
tions. 
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Fig. 1 . — Illustration of the way we measure the total size T and the separation 5 between the highest surface brightness (HSB) spots, for class a (i.e., sources 
with two HSB spots on either side of the galaxy). The contrast of these FIRST images are chosen to show various components clearly, with overlays based on the 
FIRST source catalog. The small red squares are at the locations of the sources, while the red ellipses represent the best-fit Gaussian models and enclose 95% of 
the flux. The cross shows the location of the host galaxy. In both panels north is up, east is to the left. Left: a galaxy at (a, <5) = (174.3390,61.3337) withz = 0.111. 
The image is 4' X 4'. A line that connects the two HSB spots (at the edges of the source) and the host galaxy is drawn; T is the separation between the two 
intersections of the line and the two outermost FIRST ellipses. S is the separation between the HSB spots; in this case, they fall almost along that same line. These 
two fines are shown in orange in the Figure, offset from each other for clarity. This source has = S/T = 0.93. Right: a galaxy at (a,<5) = (161.5181.38.4678) 
with z = 0.127. The image is 1.2' X 1.2'. T is the separation between the intersections of the line that passes most of the centers of the components (i.e., the 
squares) and the outermost FIRST ellipses (determined in a x^-by-eye fashion). 5 is the separation between the two HSB spots, which are the two components 
closest to the host galaxy. This source has Tj = 0.26. 



where we study the large scale clustering properties of RGs. 

Our RG catalog contains about 10,500 objects. Of 
these, 1,040 have extended morphology with roughly aligned 
lobes/jets; this is the RG sample we study in this paper. The 
details of the selection and morphology measurement are de- 
scribed in |3] Our sample is constructed to be complete in 
radio flux and optical luminosity, and is not selected against 
any particular (radio and optical) morphology, which makes it 
well-suited for investigating the FR I/II dichotomy. The min- 
imum, mean, and maximum redshift of our sample are 0.027, 
0.165, and 0.299, respectively. At z = 0.3 the resolution of 
FIRST survey (5") corresponds to a physical scale of 22 kpc. 

Finally, we cross-match all the galaxies in the NYU- 
VAGC with the DR7 of MPA/JHU-VAGCfl which pro- 
vides conti nuum-subtracted measurements of emission line 
strengths jKauffrnan n e t afl [2003 a). star-form ation rates 
jBrinchmann et a l. 2004). and stellar mass dSalim et alj 
120071), among other physical properties. These auxiliary mea- 
surements are used in ^ when we compare these properties 
among RGs and between RL and RQ galaxies. 

3. QUANTIFICATION OF RADIO SOURCE 
MORPHOLOGY 

One of our main goals is to study various properties of RGs 
and determine if their distribution is bimodal or is continu- 
ous. Our first task is therefore to define an objective measure 
(or measures) that allows us to trace the galaxy population 
smoothly from FR type I-like sources to type Il-like ones (as 
opposed to a sharp and perhaps arbitrary type I vs type II di- 

^ |http : / /www .mpa-garching .mpg . cle/SDSS/DR7/| 



vision). With the aid of such a measure, we hope to reduce 
the subjectiveness inherited in the traditional ways of classifi- 
cation, thus increasing the repeatability of our results by other 
researchers. 

Let us denote the angular separation between the highest 
surface brightness (HSB) spots on either sides of the galaxy 
as S, and the total linear size of the radio source as T . We 
follow FR and define Fs = S/T as our primary measure of the 
morphology of the radio sources. We primarily use data from 
FIRST to measure both S and T, except for very large, diffuse 
sources which become invisible at FIRST resolution. The rea- 
sons for these choices, and the details of the measurements, 
are described below. 

The quantity S depends on the angular resolution of ra- 
dio maps. It is possible that a single high surface brightness 
(SB) spot at low resolution can resolve into several spots at 
higher resolution, and the identification of the HSB spots can 
thus change. Fifteen of our RGs have been observed by a 
deep 1.4 GHz VLA survey in the SDSS Stripe 82 region (PI: 
G. Richards; see Hodge et al. 2010, in preparation), at 3 times 
higher angular resolution (FWHM ^ 1 .5") and twice the sen- 
sitivity (rms - 0.07 mJy/beam) than FIRST (FHWM 5.4", 
rms 0.15 mJy/beam). We find that the high SB spots seen in 
the deeper maps are also identified by FIRST; on average S 
measured from FIRST and from the deeper survey differ by 
less than 6%. We thus conclude that the resolution of FIRST 
is adequate to identify the highest SB spots for our sample, 
and our results should be applicable to surveys at arcsec res- 
olution for nearby RGs. More details of this comparison are 
provided in AppendixlAl 

In addition to the angular resolution, the total size T also de- 
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Fig. 2. — Similar to Fig.[T] but for class b, in which the HSB spot is coincident with thegalaxy. Left: a galaxy at (a,S) = (1 12.71 1 1,44.9336) with z = 0.072. 
The image is 1.2' X 1.2'. T is determined in the same way as the examples shown in Fig. [T] as indicated by the orange line. In this case we use the major axis of 
the central FIRST ellipse as S. This source has = 0.55. Right: a galaxy at (a, 5) = (208.6374,28.2434), z = 0.065, with r, =0.25. The image is 2.4' X 2.4'. 



pends on the surface brightness Umit of the survey. Ahhough 
NVSS is shallower than FIRST, it is more sensitive to emis- 
sion from extended sources with lo w spatial frequ ency, which 
will be absent in the FIRST maps dCondon et alj 1998). For 
most of the sources we measure T using data from FIRST; for 
41 very extended, diffuse sources (representing 3% of all the 
sources for which we have size measurement), we use NVSS 
instead. We discuss in Appendix |B] the effect of using only 
NVSS data to measure T on our results. 

For the identification of the HSB spots and the measure- 
ment of S and T, we rely on the elliptical Gaussian fits to the 
sources provided by both FIRST and NVSS (when needed) 
surveys, rather than using the atlas images, for these fitted pa- 
rameters (deconvolved size of major and minor axes, and the 
position angle) fully account for the synthesized beam proper- 
ties and are of high fidehty to the actual intensity distribution 
in the maps (iCondon et al.l[T998tlBecker"etan[T995 ). In prac- 
tice we overlay the FIRST images with ellipses based on the 
best fit parameters to guide our measurements (see Figures [T] 
and |2] for details). For each source (i.e., individual compo- 
nents of an RG), the ellipse encloses 95% of the flux, and its 
semi-major/semi-minor axes are equal to the 2a values of the 
elliptical Gaussian model. 

For galaxies associated with only one radio source (at 
FIRST resolution), or components of a multi-source RG, we 
use the fitted parameters from FIRST to determine whether a 
source is extended. A point source needs to satisfy the fol- 
lowing conditions: (1) the integrated flux-to-peak flux ratio 
fint/ fpeak < 1 -2 and (2) the deconvolved major axis < 5" (e.g., 
iBecker et a l. 1995; Kimball & Ivezic 2008). 

For every extended RG (containing in most cases at least 
two FIRST source^ whose lobes/jets are aligned to within 
~ 30° (i.e., not strongly bent as in wide/narrow angle tail ob- 
jects), we measure T and S as follows: 

' In this paper we only consider RGs associated with at least two FIRST 
sources. 



• T : (i)ln most of the cases there are two or more FIRST 
sources associated with an RG; we use the length of 
the line that passes through the FIRST source locations 
and intersects with the outermost FIRST source ellipses 
as T (Fig. [T]i. In the cases where we need to use 
NVSS for the total size measurement, we use the length 
of the line that passes through the NVSS/FIRST source 
locations and intersects with the NVSS source ellipses. 
If the size of the major axis for a NVSS source in either 
of the above cases is given as a upper limit, we remove 
the RG from our sample. 

• S: we use the peak flux of the FIRST sources to deter- 
mine the position of the HSB spots; the separation be- 
tween such spots on the two sides of the RG is S (e.g., 
Fig.[TJ. In some cases where one of the lobes is not well 
detected and modeled by FIRST, we use the FIRST im- 
ages directly. For about 1/3 of the sources the HSB spot 
coincides with the center of the galaxy (Fig. |2|i, which 
is likely due to (unresolved) jets; if the spot is extended 
and accounts for at least 10% of the total flux, we use 
the size of the spot as a measure of S. Otherwise we use 
the separation between the other high SB spots on the 
two sides of the galaxy as S. 

From visual inspection, 1,244 RGs appear to be extended 
(with apparent lobes and/or jets). We have measured S and 
T by hand for these objects using the Aladin sky atlas tool 
jBonnarel et al.ll2000l) . Based on repeated measurements for 
about fifty sources, we find our procedure is highly repeat- 
able and the resulting sizes usually agree to within 5%, except 
for very complex sources. To reduce resolution dependence 
on our classification scheme, and to remove compact sources 
(e.g., 0'Dea>1998i) . we further select a subset of 1,040 sources 
that have angular size T > 30" and physical size Tp > 40 kpc 
for the analysis presented below. Again, our RGs also satisfy 
M° 1 < -21 .27, 0.02 < z < 0.3, and /1.4 > 3 mJy. 
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TABLE 1 
Our Classification Scheme 



class 


HSBt spot(s) 


central source 


S 


number 


fraction (%) 


a 


opposite sides of RG 


not HSB 


distance between two HSB spots 


667 


64.1 


b 


center 


extended, > 0. 1/iot 


size of the HSB spot 


289 


27.8 


c 


center 


extended, < 0. l/iot 


distance between secondary HSB spots 


38 


3.6 


d 


center 


compact, > 0.3/tot 


distance between secondary HSB spots 


9 


0.9 


e 


center 


compact, < O.3/I01 


distance between secondary HSB spots 


37 


3.6 



^ "highest surface brightness" 



Motivated by the above considerations, and based on the 
ratio Tj, flux, and location of the HSB spot(s) (extended or 
point-like; coincident with the host galaxy or not), we group 
the RGs into five classes: (a) There are two HSB spots on 
opposite sides of the RG; (b) There is a single extended HSB 
spot coincident with the galaxy, with flux > 0.1/tot, where /tot 
is the total flux from all the components of the RG; (c) The 
HSB spot coincides with the galaxy, is extended, and its flux 
is < 0.1/tot; (d) The HSB spot coincides with the galaxy, is 
unresolved, and its flux is > 0.3/tot; (e) The HSB spot coin- 
cides with the galaxy, is unresolved, and its flux is < O.S/ot- 
For class b, S is the size of the HSB spot; for all other classes, 
S is the separation between the other high SB spots on both 
sides of the galaxy. 

This classification scheme is devised primarily for the ease 
of measuring r,. We do not mean to suggest there are five 
distinct types of RGs. In Table [T] we summarize this scheme, 
and list the number of sources in each of the classes. The first 
two classes (a and b) account for the majority (> 90%) of the 
extended RGs. We will focus exclusively on these two classes 
in what follows, and will examine if such an artificial classifi- 
cation corresponds to any physically different RG types, that 
is, if there is any fundamental difference between classes a 
and/7(§S|6ll. 

Does rs have any physical meaning? A HSB spot represents 
a region where the jet is highly dissipative, and thus r, can be 
thought of as an indicator of the degree of interaction between 
the jet and the environment, which depends on the strength 
and nature of the jet as well as the density of the surround- 
ing medium, among other factors. For class a sources with 

^ 1 (e.g.. Fig. [1] left), the jets can reach the edge of the 
lobes without much impediment, suggesting either a strong 
jet or tenuous ambient density, or a combination of both (e.g., 
iKawakatu et al. 2009). The opposite situation would cause a 
low rs (Fig.[ri right), or a class b object (Fig.|2]i. 

In Fig. |3] we show the distribution of for all 1,040 RGs 
(green/dotted histogram), and for the classes a (red/solid his- 
togram) and b (blue/dashed histogram) separately. Consid- 
ered together (green/dotted histogram), the distribution of 
RGs seems to be bimodal, with peaks at 0.15 and 0.60, and 
the division roughly at rj ~ 0.3. 

In the original FR scheme, the division was set to = 0.5, 
although it is not clear if FR's data showed any bimodality in 
the distribution of /-j. FR used 1.4 GHz maps with beam size 
of 23" X 23"cosec(dec) for their morphology classification. 
The difference in the division may be due to the resolution 
of the radio images (see Appendix |B]i, and/or to the sample 
used (the 3CR catalog that they used was limited to very high 
radio flux threshold and consisted of sources at much higher 
redshifts than our RGs). Nevertheless, one might be tempted 
to make a correspondence between the two distributions in 
Fig. [3] with the two FR types. 
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Fig. 3. — Distribution of r, for all the RGs (green/dotted), and for classes a 
(red/solid) and b (blue/dashed) separately. Classes c-e account for about 8% 
of the RGs in our samp le (se e Table[T). The two peaks do not correspond to 
the two FR types (e.g., j|6.U . 



TABLE 2 
Subsets of Class a 



subset 






definition 


number 


fraction (%) 


ao.9 

a<0.8 

O0.9,em 


rj > 0.8 
r, < 0.8 
r, > 0.8 
non-ao 9 


(i 
(i 
ar 
eir 


5.1 
5.1 
idL 
cla 


om > IO^-Lq (E) 
ss a objects (36. It 


86 
581 

40 
627 


13 
87 
6 
94 



Others have built upon the FR scheme and used more so- 
phisticated criteria to classify the RGs, such as the presence 
of jets or "hot spots" toward the ed ge of the lobes (e.g., OL89; 
lLeahylll993l:lGendre & Wall"2008'). In the analyses presented 
below, we will seek the best way to distinguish various pop- 
ulations of RGs; in doing so we will find ourselves defining 
several subsets of class a (aos,, fl<().8, Omaj, and oo.g.em; see Ta- 
ble |2|i. We will show that, if one sticks with a simple FR-like 
classification scheme, a division at r, « 0.8 (between what we 
call the fl<o,8 and flo.9 subsets) best separates objects of dif- 
ferent physical properties (for radio data at FIRST-resolution; 

A scheme that works even better is based both on features 
of the radio morphology (r,) and the nuclear optical emis- 
sion properties; we will refer to these as the flo.g.em (fs > 0.8 
and Lqiii > I O^Lq ) and a^aj (the rest of class a) subsets. We 
will argue in ^6.1l that the flo.g.em, flmaj, and b classes represent 
three populations of extended RGs, which roughly correspond 
to the CD, FD, and TJ types of OL89, respectively. The corre- 
spondence between these and the FR types will be discussed 
in ED 
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Fig. 4. — Distribution in the radio power-optical magnitude plane for objects in the class a (left panels) and class h (right panels). From top to bottom, the 
panels show the distribution in bins of decreasing r^. At the lower left corner of each panel we show the r, range, followed by the number of sources in that 
range. The black points (upper left panel) are class a objects with rs > 0.9, and are repeated in all other panels for comparison. There is a systematic shift of the 
locus of the RGs towards lower radio power, as r,, decreases, for class a. However, there are class a objects with lower value of r.i in the region occupied by the 
highest r, RGs. This trend is wea ker f or class b. In each panel, (green) crosses denote sources with Seyfert-like emission line ratios, while squares denote those 
with LINER-like line ratios (see Elland Table[3). 



4. BIMODALITY IN THE RADIO-OPTICAL 
LUMINOSITY PLANE? 

We start by analyzing the distribution of RGs in the radio 
power-optical magnitude plane (P-M plane for short). Owen 
and collaborators suggested that the transition from one PR 
type to the other is quite abrupt on this plane (ffl OL 89; 
lOwen & Whit"aiT99Tt IOwenl[T993l: iLedlow & Owenlfl996h . If 
this were true, and if their FR classification was solely based 
on morphological measures such as r,, we would expect RGs 
of different values of rj to occupy distinct regions in the plane. 

In Fig. ID we show the distribution of objects in classes a 
(left panels) and b (right panels). We have assumed a typi- 
cal radio spectral index of a = -0.8 (with the convention that 
/ oc v") to convert fluxes to power. From top to bottom, we 
show the distribution for RGs in bins of decreasing r,. Let us 
first focus on the class a. The black points (upper left panel) 
are RGs with > 0.9, and are shown in all other panels for 
comparison. As we will see below, many of these high r, ob- 
jects are distinct from the rest of the population in their host 
properties and environment. There is a gradual shift of the lo- 
cus of the RGs towards lower radio power, as decreases. 
Even so, the region occupied by these high-r, RGs is still 



populated by some RGs with (much) lower value of rj (e.g., 
rs ~ 0.5). Thus no particular region on this plane is inhabited 
solely by a type of object defined by some specific range of 

For class b (right panels), the trend of decreasing mean ra- 
dio power with decreasing r, is very weak, but they tend to be 
lower luminosity than the r, > 0.9 class a objects. 

These results appear to be in contradiction with those ob- 
tained by Owen and co-workers. If we were to assign an 
rs value that serves as the FR I/II divide (say = 0.9), we 
would not obtain any clean-cut separation for the RGs. Ad- 
ditional classification criteria, such as special features in the 
radio source morphology (e.g., hot spots at the edge of the 
lobes), or/and optical properties of the nuclei, may be required 
to define tw o (or maybe more) distinct populations in the P-M 
plane (see ^6.3b . 

As an independent check, we show in Fig. |5] the distribu- 
tion of FR I (blue/triangle) and II (red/circle) RGs in the P- 
M plan e using about 100 R Gs at z < 0.3 from the CONFIG 
sample (iGendre et alJl2010!) . This is a radio flux-limited sam- 
ple that was constructed by combining data from both NVSS 
and FIRST. The morphological classification was carried out 
by these authors and was based on FIRST and deeper VLA 
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Fig. 5. — Distribution in the radio power-optical magnitude plane of 
z < 0.3 FR 1 (blue/triangle) and FR II (red/circle) RGs from the flux-limited 
CONFIG sample iGendre et al. 2010), matched with the DRV NYU-VAGC. 
The morphological classification was done independently by these authors. 
There are 48 (51) type I (II) sources. The upper part (e.g., P > 10-*^ W/Hz) 
of the plane is dominated by type lis, but there is a non-negligible fraction of 
type lis at lower radio power. Mr is the (^-corrected) SDSS r-band absolute 
Petrosian magnitude. 

A-array observations. We cross-match the z < 0.3 RGs from 
this sample with the NYU-VAGC DR7 sample to obtain the 
absolute r-band magnitude. Although FR lis are on average 
more powerful than FR Is, substantial overlap between the 
two types is readily seen. This overlap still persists using a 
volume-limited subsample that consists of 34 FR Is and 28 
FR lis (sele c ted \v ith z < 0-16 and M, ^ -21. 77). In ad- 
dition, iBesj (120091) and iWin^ & Blaiitonl ( 120101) also noted 
a substantial overlap between the two FR types in the P-M 
plane. These results support our notion that in a flux-limited 
sample, there is not a sharp division among RGs in the P-M 
plane when the classification is made solely based on r,. Pos- 
sible causes of the discrepancy between our results and those 
of Owen et al. are discussed in ' 



5. HOST GALAXY PROPERTIES OF EXTENDED 
RADIO SOURCES 

We now examine the distribution of various physical prop- 
erties of the galaxies as a function of for classes a and b. 
These include the global properties of the host galaxy, derived 
quantities of the stellar populations, emission line properties, 
and properties related to the radio source and environments. 

In this section we focus on the observational results, and 
leave the interpretation t o ^ We discuss the properties of 
classes a OS.ll l and b ( ^5.2\ separately. We investigate the 
environments of the RGs in ^5.31 A comparison with RQ 
galaxies whose global properties are matched to the RGs is 
made in ^5.41 

5.1. Class a 

Throughout this paper we are concerned mainly with the 
mean behavior of samples, and thus in this section we will ex- 
amine the medians of various quantities, and their trend with 

In Fig. |6] we show the median value and its uncertaint\Fl 



The uncertainty of the median is taken as a^ /i^llN, w here u is the 
standard deviation, and N is the number of RGs (e.g., ILuptonll993) . 



for some global properties of the host galaxies as a function 
of Tj, for classes a (red) and b (blue). The small data points 
represent individual RGs, also color coded according to their 
classes. The left panels, from top to bottom, show the ab- 
solute r*^ '-band magnitude Mj.' ', rest frame color (g-r)*^ ' = 
M° ' -MJ? stellar velocity dispersion cr, Sersic index n, and 
concentration c. The panels on the right show the effective 
radius rgff, dynamical mass Mdyn, total mass density p, axis 
ratio bja for the de Vaucouleurs profile fit to the SB profile, 
and redshift. Here c = rgo /r^o is the ratio of the radii that en- 
close 90% and 50% of the Petrosian fluxes, Mdy,, = Sr^ffa^/G 
is a crude estimate of the total mass, and p = Mdyn/(47rr^ff/3) 
repres ents the (central) mass density. We follow Graham et al] 
( l2005h to estimate r^ff from «, r^o, and rgo. The r-band mea- 
surements are used for all these photometric quantities. Using 
the SMBH Mbh-ct relation (e.g., from Tremaine et al. 2002), 
we can estimate the SMBH mass Mbh; for our sample the 
range of Mbh is 2 x 10^- 1.3 x IO'^Mq, with a median of 
2.5 X IO'^Mq. 

We have also examined another measure of the galactic 
structure, the stellar mass surface density, /i = Mstar/T^j^, 
where Mjtar is the stellar mass (described below), and found 
that the trend with r, is similar to that of p. 

We focus on class a in this section. The color, Sersic in- 
dex, and concentration, and the SB profile (not shown) of the 
RGs are consistent with those of radio-quiet early type galax- 
ies. The numbers of RGs in the seven bins are (from lowest 
to highest): 118, 106, 102, 105, 101, 90, and 49. The me- 
dian redshift is almost constant for different rj bins (panel j). 
Trends or changes with are in most cases gradual and mild, 
with the exception of the highest-r, bin. Although the distri- 
butions of the properties are quite broad (c.f. the small points 
in Fig. |6l), one can identify two populations, roughly sepa- 
rated at w 0.8: sources with higher r, (hereafter the ao.9 
subset) are less luminous and massive (and thus likely har- 
boring smaller SMBH), are smaller, and have slightly higher 
density. For the rest (hereafter the fl<o,8 subset), trends with 
rs are either weak (e.g., r^fi, n, a) or absent. However, it is 
important to bear in mind that there is significant overlap in 
physical properties of RGs grouped by the r^ value. Some of 
the fl().9 objects are as m assive (cr, Mjyn) and large (rgff) as 
the fl<o.8 ones (see ^6.11 ). and, conversely, some of the latter 
population share the characteristic properties of the former. 

Four derived properties of the stellar populations (left pan- 
els), as well as emission line measurements (right panels), 
are shown in Fig. [T] All these quantities are taken from 
the MPA/JHU-VAGC. The stellar mass Ms,a, is derived from 
broadband ugriz photometry (see Salim et al. 2007 for de- 
tails). The specific star fo rmation rate (sSFR ) is der ived using 
the method presented in iBrinchmann et alj ( 120041) . The HS 
and Z)40()o indices are measures of the stellar age and star for- 
mation history (e.g., Kauffmann et al. 2003b). As in Fig. |6l 
the flo.9 and a<o.8 subsets show quite different behavior. While 
the fl<o.8 RGs have roughly the same median stellar masses, 
sSFR, and stellar ages, and do not have strong emission lines, 
the majority of ao.9 objects are quite active, but less massive, 
and have higher sSFR an d younger stellar age. For example, 
'Kauff mann et alJ (l2003al) suggested the use of [O III]A5007 
luminosity as an indicator of AGN activity; 26% of the ao.9 
RGs can be regarded as strong AGNs (e.g., logLoiii > 7). We 
do not apply any dust-extinction correction to the line lumi- 
nosities, as (1) dust content is expected to be low in the early 
type galaxies of our sample, and (2) any such correction (e.g.. 
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Fig. 6. — Median values of several physical properties of the host galaxies as a function of rj, for classes a (red) and b (blue). The errorbars show the 
uncertainty of the median values. Small points are individual RGs. Panels (a) to (j) show the absolute r^ '-band magnitude, rest frame (g — r)" ' color, velocity 
dispersion, Sersic index, concentration, effective radius, dynamical mass, mass density within r^ff, the axis ratio from the de Vaucouleurs profile, and redshift, 
respectively. Trends or changes with rj are in most cases gradual. 



TABLE 3 
Fraction of AGN-like Spectrum 







class a 








class h 




r, range 


N* 


SeyfertT 


linert 


Tj range 


Nt 


SeyfertT 


LINERT 


0.85-1.00 


46 


0.196 (0.015) 


0.043 (0.049) 


0.50-1.00 


37 


(0.005) 


(0.012) 


0.75-0.85 


85 


0.082 (0.004) 


0.024 (0.021) 


0.30-0.50 


41 


0.024 (0.005) 


(0.021) 


0.68-0.75 


91 


0.033 (0.002) 


0.022 (0.028) 


0.20-0.30 


58 


(0.002) 


0.017 (0.046) 


0.61-0.68 


95 


0.053 (0.008) 


0.011 (0.019) 


0.12-0.20 


67 


(0.003) 


(0.030) 


0.53-0.61 


88 


(0.001) 


0.011 (0.017) 


0.00-0.12 


53 


(0.008) 


0.008 (0.040) 


0.42-0.53 


93 


(0.003) 


0.022 (0.019) 










0.00-0.42 


94 


0.011 (0.004) 


0.021 (0.030) 











t Number of all RGs with good velocity dispersion and spectral line measurements. 

t The numbers in parentheses ar e for radio-quiet (RQ) galaxies whose properties are matched to the radio-loud 
(RL) galaxies in the r, bin. See j|5.4| for more details in the matching between RL and RQ galaxies. 



that based on the Balmer line ratios) is itself uncertain. 

In each panel of Fig. |4] we show as crosses the RGs 
whose spectra show signatures of Seyfert-type activity, and 
as squares those with LINER-like spectra. We follow 
iKauffmann et alj (I2003al) to classify the spectra according to 
the [O III1A500 7/H/3 and [N II1A6583/Ha hne ratios on the 
BPT diagram ("Baldwin et aP I198U IVeilleux & Osterbrockl 
[T987h . The class a sources with higher values of rj contain 
larger fractions of active nuclei. For example, 6/21 « 29% 
of the r, > 0.9 sources have a Seyfert-like spectrurrFT, and 
2/21 « 10% can be classified as LINERs. In Table [5] we 
record the fraction of RGs with Seyfert- and LINER-like nu- 
clei as a function of r,. Note that the r, binning in the Table is 

" Although there are 24 RGs with r, > 0.9 (Fig.|4), only 21 have spectra 
of good enough quality for line measurements. 



the same as that in Figures|6]-[8] chosen such that each bin con- 
tains roughly equal number of RGs (except for the highest 
values), and is thus different from that shown in Fig.lD While 
the fraction of class a RGs with LINER-like spectra is roughly 
the same for all bins (at about few percent), the fraction of 
Seyfert-like spectrum is a strong function of r^. About 20% 
of the objects in the aq 9 subset have strong emission lines that 
are characteristic of active nuclei. 

As an indicator of the accretion rate onto the SMBH, we 
show in Fig.[8](panel g) the [O III] line Eddington ratio, which 
is the [O III] A5007 luminosity divided by the Eddington lumi- 
nosity LEdd (which is estimated using the Mbh-c relation). 
The median accretion rate is quite close to zero for most of 
the rs bins, but is higher for the oo.g systems. In panel (f), we 
show the analogous radio Eddington ratio (Pi,4/LEdd) at 1.4 
GHz, which exhibits the same trend as Loiii/^-Edd- Note that 
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Fig. 7. — Median values of several derived properties of the host galaxy stellar populations, as well as emission line measurements, as a function of r,, 
for classes a (red) and h (blue). Panels (a) to (h) show the stellar mass, specific star formation rate, US index, 4000A break strength D4000. [O III]A5007 line 
luminosity (no extinction correction applied), [O III] equivalent width (EW), Ha line luminosity, and Ha EW, respectively. The errorbars show the uncertainty 
of the median values. Small points are individual RGs. 



we do not integrate over the radio spectrum (e.g., from 0. 1 to 
10 GHz) to obtain the total radio power; but as we are mainly 
interested in the trend with rj, this should not be a problem 
under the assumption that the spectral shape is not a function 
of rs- Assuming a mean spectral index of a = -0.8 for all RGs 
and integrating over the frequency range = 0. 1 - 10 GHz, we 
would need to scale the 1.4 GHz Eddington ratio by a factor 
of 74 to obtain the radio Eddington ratio. 

The other panels in Fig. [8] are related to the environments 
and other radio properties of the RGs. As a simple, statistical 
measure of the number of neighbors, for each RG we count 
the number of luminous galaxies (M* - 3 .5 < Mj?' ' < M* + 1 .5) 
within 1 Mpc and 0.5 Mpc in the SDSS photometric cata- 
log, assuming they are at the redshift of the RG, and sub- 
tract the expected number of galaxies (from the global galaxy 
counts) in the same apparent magnitude range. These neigh- 
bor counts are denoted as Ei and E0.5, respectively, and are 
shown in panels (a) and (b). These two scales are chosen to 
reflect the scales of the host group/cluster and of any possible 
local structure within the group/cluster. On average, class a 
RGs live in dense environments (i.e., in excess with respect 
to the global background). The median number of neighbors 
decreases with increasing r,. 

Panel (c) shows the mean 1 .4 GHz radio power as a func- 
tion of r,. Here objects in the flo.9 subset have the highest 
luminosity, while a<().8 RGs have very similar median radio 
power. Panel (d) shows the distribution of physical sizes of 
the RGs. Not surprisingly, the RGs with highest r, are also 
largest (median size exceeding 300 kpc), since it is mainly 
the lobes that contribute to the low-frequency radio fluxes. 

Panel (e) is another measure of the radio source morphol- 



ogy: the central-to-total flux ratio r/. We estimate the central 
flux by summing fluxes from all components within O.IST 
from the center. For class a, the ry ratio is on average small 
(< 0.1). Note that this must be partially a selection effect, for 
if a strong source is present at the center, it will likely be the 
HSB spot and the RG would be put in other classes (b-e). 

Finally, panel (h) shows the logarithm of the radio- 
to-optical luminosity ratio (calculated simply as log f 1.4 + 
0.4M° '), which is also almost constant at rj < 0.8, and rises 
sharply at highest rs bin; this enhancement for aq 9 is due to 
the combined effect of it having the lowest median optical lu- 
minosity and highest median radio power among the class a 
objects. 

All the results from Figures |6]-[8] indicate that there appear 
to be two populations of RGs among class a: galaxies in the 
ao.9 subset have lowest stellar mass, but have the highest star 
formation rate, AGN activity, and radio power, and live in rel- 
atively poorer environments, while the fl<().8 RGs, especially 
those with r, « 0.4-0.8, share very similar properties, such 
as the luminosity, size, dynamical mass, density, indicators of 
recent star formation history (e.g., D4000 and H6 indices), and 
neighbor counts within 0.5 and 1 Mpc. We should emphasize 
that the division in rs of the two populations is not sharp; we 
choose to distinguish the two subsets at = 0.8 mainly for 
simplicity. 

Can we identify these two populations with the two FR 
types? The general properties of these two populations do 
conform roughly to the characteristics of the two FR types 
found in the literature (|T]i. Before making a direct correspon- 
dence (see ^6. Il l, however, it is important to investigate where 
the class b objects fit in the context of the FR I/II dichotomy. 
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Fig. 8. — Median values of environmental and radio properties of the host galaxies as a function of rj, for classes a (red) and h (blue). Panels (a) to (h) show the 
excess number of neighbors over the mean background within 1 Mpc and 0.5 Mpc (Si and So. 5), the 1.4 GHz radio power, the linear physical size of the radio 
sources, the central-to-total radio flux ratio, the radio and [O III]A5007 Eddington ratios, and the radio-to-optical luminosity ratio, respectively. The eiTorbars 
show the uncertainty of the median values. Small points are individual RGs. 



5.2. Class b 

The r, binning in Figures |6]-[8] for the class b (blue points) 
is the same as that shown in Fig.|4](right panels); the numbers 
of RGs in the five bins are (from lowest to highest): 57, 78, 
68, 44, and 42. The median redshift slightly increases with r, 
(Fig. in panel j; see Appendix IaIi. 

Broadly speaking, the class b RGs are typical early type 
galaxies in their optical properties (based on e.g., the color, 
concentration, and Sersic index). One gets an impression 
from Figures |6]-[8] that there is not much variation with rj for 
galaxies in class b. For example, reff, Mdyn, p, /i, and 1.4 
GHz Eddington ratio all have median values that do not vary 
much with r,. The most notable systematic variations are that 
(g-r)" ' color becomes redder as r, increases, a positive cor- 
relation between r/ and r^, and a negative correlation between 
the total size of the radio source and r^. 

There is considerable overlap in the distributions of various 
physical properties between the fl<o.8 and class b RGs. How- 
ever, the class b objects have slightly higher median a, c, p, 
p,, smaller reft, lower emission line strength, neighbor counts, 
radio power and Eddington ratio, and radio-to-optical lumi- 
nosity ratio. That is, they are more compact and are more qui- 
escent in terms of their nuclear and radio activity. Of the 597 
class a objects for which classification of spectral properties 
using the BPT diagram is possible, 37 (= 6.2%) are Seyferts 
or LlNERs. For class b, the fraction is only 4/256 = 1 .6%. 

In Table |4] we list the basic properties of the RGs in our 
sample. 



5.3. Environments of Extended Radio Galaxies 



From Fig. |8] we see that the £3<o.8 subset has the highest 
neighbor counts (median S0.5 ~ 9.9 ±0.6), while those of 
class b and 00.9 are close (S0.5 ~ 7.4 ±0.6 and 6.7 ± 1.1, 
respectively). For reference, for the most luminous galaxies 
(MJ.'-' = -22.77 to -24.27, irrespective of radio properties), 
E0.5 ~ 10.8 ±1.0. 

To better understand the environments of these RGs, we 
further examine their association with galaxy clusters, and 
calculate their clustering properties. First, we cross-match 
our RG sample with the clu sters found by the maxBCG al- 
gorithm (iKoester et al.ll2007h . 704 (= 75%) of our class a and 
b RGs lie in the footprint and the redshift range (z = 0.1-0. 3) 
of the maxBCG survey. We focus on clusters more mas- 
sive than K, 1O''*M0, above which the cluster catalog is about 
~ 80% complete. The cluster mass and the virial radius are 
estimated using the w eak lensing-calibrat ed mass-observable 
scaHng relation from iReves et all (120081) . Only 13%, 24%, 
and 14% of the RGs in the ao.9, a<o.8, and b subsets are as- 
sociated (i.e., within one virial radius) with these optically- 
selected cluster£3- However, > 80% of those RGs that are 
in clusters lie within 20% of the virial radius from the cluster 
center, consistent with previous find ings that RGs are cen- 
trally concentrate d in clusters (e.g., iLedlow & OwenI [19951; 
iLin & Mohiil2007l) . 

The RGs are mainly hosted in dark matter halos more 

The class a and b RGs only account for < 10% of our full RG sample. 
The narrow-angle tail objects, whose morphology is believed to be due to the 
relative motion of the galaxies and the intracluster medium, together with the 
compact sources (which account for the majority of the sources), are excluded 
in the current sample. The fraction of all RGs that are cluster members is 
therefore higher than we have estimated here. 
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Fig. 9. — Redshift space RL-RQ cross correlation functions for tlie a^ g 
(black triangles), a<o.8 (red circles), and b (blue squares) subsamples. These 
are calculated using galaxies that form a volume-limited sample (z < 0.16, 
MJ" < -21.77; for RGs an additional criterion is logP > 23.31). At scales 
larger than K, 200/;-' kpc, both a<o.8 ™d class b have similar clustering 
strength, suggesting that they are hosted by halos of similar mass. The ao.9 
subset might be hosted by halos of lower mass, although a larger sample is 
needed to test this possibility. Some a^ g RGs have similar optica l pro perties 
and likely also similar environment as the other two subsets (see 86. It . 



massive than « 2 x lO'-^M© (iMandelbaum et all l2009h . so 
the small fraction of RGs associated with maxBCG clus- 
ters suggests that the majority of our RGs must be associ- 
ated with clusters or groups of mass ss 2- 10 x lO'^'M©, a 
range in which the maxBCG catalog is highly incomplete 
(iKoester et al.ll2007l) . We then resort to the clustering proper- 
ties of the RGs, which provides some insight into the relative 
mass scales of halos that host these subsets of RGs. We cal- 
culate the cross correlation functions between the RG subsets 
and the general galaxy population, constructed as a volume- 
Umited sample of 73,202 galaxies from the NYU-VAGC DR6 
(z < 0.16, M," ' < -21.77). We select the RGs to satisfy the 
same redshift and magnitude cuts as well as a lower limit in 
radio power logPi 4 > 23.31, which results in 41, 211, and 
123 objects for oo.g, fl<o.8, and b subsets, respectively. The 
redshift space cross correlation functions are shown in Fig.|9] 
for flo.9 (black triangles), fl<o.8 (red circles), and class b (blue 
squares). Although all three subsets have similar clustering 
strengths at scales > 200/1"' kpc, there is a slight hint of lower 
clustering amplitude for the aq 9 objects; if this is confirmed 
with larger RG samples, this implies the host halos of the ao.9 
subset are on average less massive than the hosts of the other 
subsets. 

Combining these results, we see that the distributions of 
host halo mass for the a<o.8 and class b subsamples are similar 
at group mass scale, but that of the former must have a higher 
tail towards clusters (> lO'^M©). For RGs in groups, the fl<o.8 
objects must be more centrally concentrated than the class b. 

5.4. Radio-Loud vs Radio-Quiet 

So far we have made comparisons among different subsets 
of RGs. It is important to place them in the context of general 
massive galaxy populations. Furthermore, some of the cor- 
relations between the physical properties and r, may be due 
to other fundamental correlations of the early type galaxies 



that have no physical connection to the radio source. For ex- 
ample, given that optical luminosity and e ffective radius ar e 
tightly correlated in elliptical galaxies (e.g..'S hen et al.ll2b03h . 
an anti-correlation between Mr and r, would imply a simi- 
lar anti-correlation between r^ff and r, (Fig. |6ll. To take out 
such an effect, for every RG we find up to ten RQ galaxies 
that have very similar redshift and absolute magnitudes in the 
^o.i^ r*' and /*' ' bands. For all RGs in a given bin, we 
calculate the average value of the physical property in ques- 
tion from all the matched RQ galaxies and subtract that value 
off from the mean obtained for the RGs (or take the ratio of 
the two, depending on the nature of the properties). Note that 
some AGNs (Seyferts and LlNERs) may be included in the 
RQ sample, as we do not distinguish truly quiescent galaxies 
from those that can be regarded as AGNs based on the optical 
emission line diagnostics. 

Figures [To] and [TT] show the results as a function of r,. The 
comparisons presented in these two Figures are between each 
class and its RQ counterparts, not between the two classes. 
For most of the properties shown in Fig. [TO] (environments 
and global galaxy properties), RL galaxies have higher mean 
values than do the RQ ones (e.g., more neighbors, higher mass 
and velocity dispersion, larger effective radius). The only ex- 
ception is the central mass density/stellar surface density, for 
which the fl<o.8 subset appears to be less dense than its RQ 
counterparts, mainly due to their larger scale lengths. Fig.fTTI 
shows comparisons for some derived quantities of the stellar 
population and the line measurements. The fl<o.8 and class b 
RGs on average have slightly higher stellar mass, and older 
stellar age (AHJ = H(5rl-H(5rq < 0), than their respective 
RQ counterparts. 

Among the class a objects, galaxies in the aq 9 subset are 
closest to their RQ matches in terms of neighbor counts, reff, 
Mdyn, and /.t, but show dramatically stronger emission lines 
and sSFR. In contrast, emission lines are weaker for class 
b RGs than their RQ counterparts. Class b objects also live 
in environments that are closer to their RQ counterparts than 
those of class a, and their sizes (^eff) are also more similar. 

Let us examine the spectral properties of the RGs in more 
detail. In Table [3] in parentheses, we list the fraction of 
the matched RQ galaxies that show Seyfert- and LlNER-like 
spectrum. This is useful for evaluating whether the active nu- 
cleus fraction is elevated in a given bin. For example, for 
class a with > 0.85, 19.6% of RGs have Seyfert nuclei, 
while only 1.5% of the matched RQ galaxies exhibit the same 
level of activity. In the same r, bin, 4.3% (4.9%) of the RL 
(RQ) galaxies have LINER nuclei. The LINER fraction for 
class a is roughly independent of rs, and is fully consistent 
with what is found in the RQ populations, but the Seyfert ac- 
tivity is enhanced by about a factor of ten in RGs with rj > 0.6 
(see also llvezic et al.l2002l) . There is some suggestion that the 
nuclei of class b galaxies are actually more quiescent than the 
matched RQ galaxies. For example, there are 2659 RQ galax- 
ies matched to the class b subsample, out of which 84 are 
LlNERs. We thus expect to find about 8 LlNERs in the 256 
class b RGs, but only detect 3, w hich is inconsis tent with the 
Poisson expectation at 99% level (lGehrelslll986l) . 

The differences in the optical spectral features are further il- 
lustrated in Fig. [12] which shows stacked SDSS spectra of var- 
ious subsamples of RGs. To make a fair comparison among 
the three subsets (00.9, fl<o.8, b), as well as between RQ and 
RL galaxies, we select galaxies in three SMBH mass (stellar 
velocity dispersion) bins, and limit the ranges of dynamical 
mass and stellar mass surface density to be within 50% of 
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Fig. 10. — Differences between RL and RQ galaxies. Red and blue points denote classes a and b, respectively. For every RG we find up to ten RQ galaxies 
that have very similar M° KM°-', and M°- ' magnitudes and redshift; after calculating the mean value of a given physical quantity from all RGs in a given bin, 
we then take the ratio (or the difference if the quantity is logarithmic) between the mean value of the RGs and that obtained from all the RQ galaxies that are 
matched to the RGs in question. The left panels, from top to bottom, show the ratio in neighbor counts within 1 Mpc and 0.5 Mpc, the velocity dispersion, Sersic 
index, and concentration. The right panels, from top to bottom, show the ratio/difference in effective radius, log dynamical mass, log mass density, contribution 
of the de Vaucouleurs profile to the SB profile (normalized to unity), and the axis ratio. The errorbars denote the uncertainties in the mean value. 



the flo.9 locu£3- These properties are chosen to select RGs 
of similar central engine, fuel supply, and structure. There are 
(19,46,23) RGs in the lowest mass bin for the fl<o.8, and b 
subsets, respectively. In the intermediate and high mass bins, 
the numbers of RGs ai-e (1 1,99,63) and (1 1,71,40). In the Fig- 
ure, the three columns correspond to the three mass bins (in 
increasing mass order from left to right); three pairs of panels 
are shown in each column (from top to bottom: oo.g, fl<o.8, 
b). The pair consists of the mean spectra of an RG subsam- 
ple and the RQ galaxies of similar Mbh, M^y^, and /i (upper 
panel) and the differences between the RL and RQ spectra 
(lower panel). For each Mbh bin, the RQ galaxy mean spec- 
trum is an average over 80 randomly selected galaxies, which 
are chosen irrespective of their spectral properties, and thus 
may contain some AGNs (Seyferts and LlNERs). 

A few points are worth noting from the RL-RQ spectral 
difference panels. (1) Only 00,9 RGs have statistically signifi- 
cantly stronger emission lines than their RQ counterparts. (2) 
The emission line (particularly [O III]A5007 and Ha) strength 
decreases as Mbh increases for ao.9 objects. (3) In the lowest 
mass bin, a<o.8 and class b RGs are redder than their RQ coun- 
terparts (based on the difference spectrum); this is not seen in 
the other binsQ In addition, the Ha line is much weaker than 
that in the RQ galaxies. (4) In each mass bin, fl<o.8 and class b 

For computing the mean spectra, the RQ galaxies are selected to match 
a, Mjyn, and /i of the RGs, without any restriction on redshift and restframe 
g" ', r , and ' magnitudes, as was done eariier in this section. 

Recall that in Figures [Tol and rm the optical colors [(g — r)" ', (r-f)" '] 
are used to match the RL and RQ galaxies, and thus we do not examine color 
differences. 



have, to first order, similar difference spectra, indicating their 
spectra are close to each other. 

6. DISCUSSION 

In this section we build upon the observational results pre- 
sented in the previous two sections to investigate some in- 
triguing questions related to the physical nature of RGs: What 
is the correspondence between the FR types and the classes a 
and b? Do class b or a<o.8 objects represent evolutionary 
sequences? What is the physical origin of various morpholo- 
gies ? 

6.1. Three Types of RGs? 

46 out of 85 flo.9 objects have no detectable [O III]A5007 
emission line (i.e., the signal-to-noise ratio of the line is less 
than three). How do they differ from those objects with emis- 
sion lines? About 75% of 00.9 RGs with [O III] line luminosity 

> IO^Lq (roughly corresponding to [O III] Eddington ratio of 

> 10"^) show clear hot spots at the edge of the lobes (giving 
the impression of a bullet shot into a tenuous medium), while 
about 2/3 of the 00.9 objects without emission lines have HSB 
spots that show less contrast with the lobes, or have lobes that 
are not well-aligned. It is plausible that the mechanism that 
creates the emission lines is physically related to the process 
responsible for the generation of hot spots. 

The distributions of many physical properties for the ao.9 
RGs with and without emission lines are often offset from 
each other (albeit with substantial overlap). It is the oq 9 RGs 
with emission lines (hereafter ao.9,em objects; Table |2|i that 
make this subset stand out from the rest of RGs. On the other 
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hand, the flo.9 objects without emission lines share very simi- 
lar properties with the aq s subset; the median values of most 
physical properties are within Icr of each other We may re- 
gard them as the same population as the fl<().8 subsample, and 
refer to the combined population, which accounts for the ma- 
jority of class a objects, as Omaj (see TableO. 

We emphasize that the ao.g.em objects are still massive 
galaxies (M,- < M*), and that the distinction of this popula- 
tion from the rest of RGs is far less dramatic than the red- 
blue galaxy bimodality of the general galaxy population (e.g., 
Haldry et al. 2004). 

Since our proposed classification scheme combines both 
radio morphology and nuclear emission line strength, while 
previous ones usually rely on one or the other of these crite- 
ria (see a perfect correspondence between the two is not 
expected. We also emphasize that our scheme is more quan- 
titative, objective, and reproducible, than are either the FR or 
OL89 classifications. 

Broadly speaking, we can identify the flo.9,em objects with 
HE RGs, and the rest in our sample with LE RGs, with the 
caveat that some of the ciniaj objects do have strong emission 
lines (c.f. Table O. Regarding the classification scheme of 
OL89, the flo.9.em subset corrcsponds to the CD type (mainly 
because of the high occurrence of hot spots of the former), the 
flmaj subset is consistent with the FD type, and class b coin- 
cides with the TJ type. Although the average optical proper- 
ties of class b and fl<o,8 (or Omaj) RGs are similar, we regard 
them as distinct populations, mainly based on the differences 
in their environment, radio, and nuclear activity. 

By definition, the flo.9,em objects have > 0.8, and are thus 
associated with FR lis according to the original FR definition. 
A direct correspondence between the FR types and our flmaj 
and b subsets is not possible, however, given that the distri- 
butions for both class b and flmaj objects are quite broad (e.g.. 
Figures |6]-[8ll. If we have to adopt a dichotomy classification 
scheme, as advocated by FR and followed by many others, 
then we may call the class b plus fl<().8 objects type I, and the 
ao.9 objects type 11. However, referring to these two groups as 
"FR 1" and "FR 11" would be mis-leading, as the original FR 
definition is solely based on r,. 

To summarize, based on the properties of the host galax- 
ies and radio emission, we suggest there are three groups of 
RGs in our sample: ao.9,em, Qmaj, and class b. There is no sin- 
gle physical property that can be used to cleanly separate one 
group from the others. For example, the distribution of the 
flmaj group is quite broad (c.f. Fig. [Sj, almost encompassing 
that of the ao.g.em objects at the high r, end. A simple mor- 
phological measure such as rs is thus only of limited use for 
classifying extended RGs. 

In addition to r^, we have also explored the use of radio 
power Pi 4 in the classification scheme. Although the ^ 5% 
of the RGs with highest radio power are reasonably separated 
from the rest in plots like Figures |6H8]( with abscissa replaced 
with f 1.4), suggesting that selecting via Pi .4 can in principle 
producing a subsample similar to ao.9, we decide to stick with 
the morphological parameter r^, as understanding the origin 
of differences in the radio morphology is one of the main ob- 
jectives of this paper (see ^6.4t . 

We conclude by estimating the abundances of the three 
morphological groups, using the same volume-limited sample 
as in E3](z < 0.16 andM° ' < -21.77; for RGs an additional 
requirement is logPi 4 > 23.31). The abundances relative 
to all galaxies are (flo.9,em,flmaj,^) = (0.034%,0.41%,0.20%). 
Among the RGs (irrespective of morphology/extendedness). 
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Fig. 13. — Illustration of a possible evolutionary sequence for class h RGs 
(cyan open circles). These objects are selected to have MJ? Mjyn, /i, and 
So. 5 similar to the median value of class b RGs in the rs = 0.2-0.3 bin. The 
blue solid points and the errorbars are the median value and its uncertainty for 
all class b objects. It is plausible that during the life of a class b source, both 
f-j and rf decrease with age, while Tp increases, as implied by the average 
trends of the cyan points. The magenta squares are the fl<o,8 objects with 
similar mass, structure, and neighbor counts as the cyan points. They do 
not follow the trends of the cyan points, suggesting that for this particular 
choice of mass and environment, these class a objects do not represent the 
later phases of evolution for class b objects. 

the fraction of these types are 0.9%, 11.1%, and 5.5%, re- 
spectively (see also Table|2]). 

6.2. Evolutionary Sequences? 

The small dependence of the majority of physical proper- 
ties we have examined on for the class b objects prompts 
the question: are they the same RGs viewed at different stages 
of evolution? One could imagine that a young RG starts with 
large and rf (and smallest total size Tp); as the jets/lobes 
advance, and rj both decrease, while Tp grows. Such a 
trend seems to be present in Fig. [8] It is important to realize 
that, however, many such evolutionary sequences (of differ- 
ent combinations of the host galaxy, central engine, environ- 
ments, etc) are probably simultaneously present in our sam- 
ple, and therefore the median behavior of the class (as seen in 
Figures |6H8]) may not reflect any one sequence. To single out 
an evolutionary sequence, one should therefore only consider 
RGs of very similar properties (at least those properties that 
will not change over radio lobe time scales), such as mass, 
structure, and neighbor counts. 

We test this idea in Fig. [13] The cyan open points are a 
subset of galaxies in class b selected to have M, , Mdyn, M, 
and S().5 similar to the median value of the = 0.2-0.3 bin. 
If they can be regarded as an evolutionary sequence param- 
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Fig. 14. — Distribution of extended RGs in tlie P-M plane. RGs without 
appreciable [O III]A5007 emission are shown as small green points. Those 
that do are represented as larger points whose color reflects the [O III]A5007 
Eddington ratio (OER), as indicated on the color bar on the right hand side. 
flo.9 objects without and with OER measurements are shown as triangles in 
magenta and in orange of various degrees of intensity, respectively. 



eterized by decreasing r,, we would expect r/ to decrease 
while Tp increases, which is in rough agreement with the ob- 
served trends (of the open points), although we caution the 
trends may be somewhat driven by the few objects at larger 
r,PI Even though their size becomes bigger, their radio power 
stays about the same, probably due to the fast dissipation of 
energy in the jets; the outer regions do not contribute much to 
the luminosity^ 

We saw in ^5.2l that there is significant overlap in the distri- 
butions of various properties for class b and a<o.8 RGs. It is 
possible that some of the fl<o.8 RGs represent the later phases 
of evolution of class b objects. In Fig.[T3]we show as magenta 
open squares the fl<o.8 objects with the same ranges of Mj? 
Mdyn, and S0.5 as the class b objects (cyan points). The 
most notable trend with r, is Tp. If these subsets of the two 
classes were related, rj needs to increase as the sources age. 
However, the higher typical radio power of the fl<o.8 objects 
make such an evolutionary scenario implausible. 

6.3. The P-M Plane Revisited 

In ^ we noted that RGs in our sample are not separated 
into two groups in the P-M plane via a simple division in 
Tj. With the correspondence between our three subsets and 
the three morp hological groups identified by Owen and co- 
workers ( ^6. lb , could we better reconcile their results with 
ours? 

In Fig. [14] the small green dots show all the extended RGs 
in our sample, and triangles are the aq 9 subset. The color bar 
in this Figure represents the [O III] A5007 line Eddington ratio 

Another caveat is that the surface brightness of the class b objects may 
decrease rapidly outwards, producing correlations between r,, Tp and ry sim- 
ilar to that due to an evolutionary sequence for sources close to the detection 
limit. We thank Philip Best for pointing this out. A larger sample, with care- 
ful selection criteria, would be needed to assess the contamination due to this 
effect. 



(OER). There is a better, although still not complete, separa- 
tion in the P-M plane of the flo.9,em RGs with relatively higher 
level of nuclear activity (e.g., OER > 10"^, roughly corre- 
sponding to Loiii > lO^L©) from the rest of the population. 

The significant overlaps among the different classes of RGs 
in the parameters we have surveyed (e.g.. Pi 4, MJ? ', mass, 
structure, environment; Figures |3-[S) imply that no simple 
combination of non-radio observables can be used to deter- 
mine/predict the morphology of the RGs (which may suggest 
there are many physical processes that determine the radio 
morphology). It is also possible that the overlap in properties 
is due to a mixture of RGs at different stages in their evolu- 
tion, as we argued above. 

We suspect the discrepancy between our finding and that 
of Owen et al. - who stated that the FR I and II RGs could 
be separated "cleanly" on the P-M plane - is due to sam- 
ple construction. Substantial overlap between the two types 
is apparent iising e ither our sample or the subsample from 
iGendre et al.l ( 2010 , see Fig. [5]l, both are r adio flux-limited 
(see also lBesi3l2009l:IWing &Blantonll2010l) . However, sam- 
ple selection was quite heterogeneo us in some of the earlier 
works (OL89: lOwen & WhitdlloMl IOwenlll993h . where the 
main criterion for inclusion was to sample the P-M plane 
as much as possible. In this sense, FR II objects are over- 
represented. With samples assembled under better defined 
criteria (e.g., limited by flux and redshift, restricted to central 
parts of clusters), some overlap between the two types were 
seen in Ledlow & Owen (1996). 

6.4. Origin of Different Morphologies? 

The existence of objects whose two lobes exhibit clearly 
differ ent FR mor phologies from each other ("HYMORS" ob- 
jects, iGopal-Kris hna & Wiita 2000) argues that the immedi- 
ate surrounding of the host galaxies must play some role in 
shaping the morphology. In addition, analytic models of ra- 
dio sources suggest that a key quantity in determining the 
large scale morphology is the jet pow er-to-ambi ent density 
ratio (Ljl Pn) of the h ost galaxies (e.g., iKaiser &" Best 2007t 
'Kawakatu e t al.l l2009i and references therein). Here pa is 
measured at the core radius of the host galaxy. Jets from hosts 
with low values of Lj/ pa are more prone to the development 
of turbulence and become subsonic, resulting in plume-like 
morphology beyond the deceleration point, while jets from 
systems with high Lj/ pa ratio are strong enough to remain 
supersonic, leading to the hot spots at the edge of the lobes. 
Since the jet mechanical power is not directly observable, we 
assume it is proportional to P^ J (e.g., Cavagnolo et al. 20ldl 
and references therein). We also assume p can be regarded 
as a faithful proxy for the local interstellar medium density 
in the host galaxies. The median mechanical power-to-total 
mass density ratio (cx Pf-J j p) for oo.g.em, Qmaj, and b is roughly 
3.2 : 2 : 1, a trend in qualitative agreement with the models. 
In fact, the median values of p for the three subsamples vary 
only by about 35%, and it is mainly the difference in the radio 
power that drives the Pi'H p ratio Pi 

We suspect that the accretion rate onto the SMBH is more 
important than the structure of the galax y in determining the 
differ ent radio morphologies (see also Bau m et anil995l : iHol 
120081; f or discussi ons on the modes of a c cretion, see e.g., 
' Best et al., 2005bt iHardcastle et"!!] 120071: iKauffmann eTall 
120081) . At high accretion rates, the primary accretion flow 
is likely a geometrically thin, optically thick disk, which may 

The results are qualitatively the same if we use the P." ? //i ratio instead. 
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launch jets that are very well-collimated over hundreds of kpc 
scale. As the accretion rate decreases, the thin disk moves 
away from the SMBH, and the inner region is o ccupied by 
a radiatively ineffici ent accretion flow (RIAF; see lEsin et ali 
ll997tlNaravanl2 005'). with diminished emission line luminos- 
ity. Such an accretio n flow is known to create outflows (e.g., 
iNaravan & ^119951: Blandford & Begelman 1999), and it is 
likely that jets so created will not be well-collimated (e.g., be- 
yond hundreds of kpc), or are collimated initially, but suffer 
entrainment and deceleration very early on due to their lower 
intrinsic power. 

If this picture is correct, we may understand the three sub- 
sets as follows. First, ao.9.em objects have highest SMBH 
accretion rates (e.g., OER > 10"^), are powered by classi- 
cal thin accretion disks with strong, well-collimated jets that 
can produce strong hot spots, and are usually associated with 
lower mass galaxies living in less dense environments (with 
respect to the other subsamples considered here). The Omaj 
RGs are massive, found in dense environments, and their cen- 
tral engines are likely fed by lower accretion rates (e.g., OER 
< 10" ), probably in a RIAF. Finally, those galaxies with low 
accretion rates and with relatively low Lj/pa ratio will likely 
show jet-dominated morphology, making them class b RGs. 

To some degree the accretion rate correlates with r^; vari- 
ations in the SMBH spin, magnetic fields, the structure of 
the galaxy, and the density of the intergalactic/intracluster 
medium, however, may all cause spreads in rj at a given ac- 
cretion rate. It is possible that uncertainties in Mbh (inferred 
from the Mbh-c relation) also smear the correlation. About 
12% (9%) of the fl<o.8 (class b) RGs have OER > 10-^ while 
these high accretion rate class b objects may represent ear- 
lier phases in evolution of CDs, the active class a RGs with 
Tj < 0.8 may be manifestations of the variations in at a given 
accretion rate mentioned above. 

The difference in the accretion rate (and in turn the jet- 
launching mechanisms) may also explain the spectral prop- 
erties of the host galaxies (see Figures [7] and [12] |. For jets 
created by a thin disk, the "zone of influence" within the host 
galaxy is quite small (e.g., the jets may only punch two small 
"holes" in the galaxy), and thus any feedback due to the jets 
cannot suppress efficiently the star formation activity that may 
be linked to the onset of the AGN activity. On the other hand, 
if the jets launched by a RIAF are not well-collimated, they 
may influence a much larger volume of the host galaxy and 
thus terminate star formation more easily. The most efficient 
feedback mode (for the host galaxy itself) may be a combina- 
tion of a RIAF and a dense interstellar medium (or immediate 
surrounding of the host galaxy), which slows down the jets 
quickly and creates the class b morphology. 

7. CONCLUSION 

Extended radio galaxies (RGs) have been classified based 
on their radio morphology or nuclear emission line activity. 
In this paper we have proposed a hybrid classification scheme 
that combines both features, and presented a comprehensive 
study of the host galaxy properties of RGs. Our main ob- 
jectives are to detect and define distinct populations of RGs, 
to understand the traditional Fanaroff-Riley (FR) type 1/11 di- 
chotomy in the context of our new scheme, and to unravel the 
origin of different radio morphologies. Our RG sample con- 
sists of 1,040 objects selected with 1.4 GHz radio flux density 
/i.4 > 3 mJy, r'^-^-band absolute mag nitude M° ' < -21.27 
(i.e., more luminous than M*, the characteristic magnitude 
of the galaxy luminosity function), radio angular diameter 



T > 30", physical size Tp > 40 kpc, and at 0.02 < z < 0.3. 
All of the RGs in our sample appear to be massive early type 
galaxies. 

We use the [O III]A5007 line luminosity as an indicator of 
the nuclear emission strength, and use a parameter = S/T 
as a continuous parameterization of the RG radio morphol- 
ogy. Here T is the total size of the radio sources, and S is 
the separation between the highest surface brightness (HSB) 
spots on either sides of the galaxies. Roughly 60% of our ob- 
jects show HSB spots on both sides of the host galaxy; we 
refer to these as class a RGs. About 30% of the sources ap- 
pear to have prominent jets, with HSB spot coincident with 
the host galaxy. We call this population class b (see Figures[T] 
and|2l Table[I]i. 

Our main results are as follows: 

• The distribution of is bimodal (Fig. [3]), although we 
argue that t he tw o peaks do not correspond to the two 
FR types (SM. 

• Among the class a objects, a small population with high 
values of (> 0.8) andhigh [O III] A5007 line luminos- 
ity (Loiii > IO^Lq) seems to be distinguished from the 
rest, in the sense that on average they are hosted by 
lower mass galaxies, live in relatively sparse environ- 
ments, and have higher accretion rates onto the central 
supermassive black hole (SMBH), as manifested by the 
[O III]A5007 line Eddington ratio (EH Figures IgHS]). 
We refer to these RGs as the oo.g.em subset, and the rest, 
the majority of class a, as the Omaj subset (Table|2|. The 
distribution of r, for the a^aj objects is quite broad, en- 
compassing the range occupied by the flo.g.em RGs at the 
high end. A simple morphological measure such as 

is thus only of limited use for classifying extended 
RGs. 

• The average properties of flmaj and class b RGs, such 
as the optical luminosity, stellar mass, 4000A break 
strength, and velocity dispersion, differ by 20% or less 
from one another. However, because of the differences 
in the environments (e.g., characterized by the neigh- 
bor counts within 0.5 Mpc; Fig. [8] panel b) and the 
(nuclear) emission line properties (Fig. [8] panels g & 
h; Table [3 ]t, we regard them as distinct populations 

(§£2E1. 

• Among the three subsamples (flo.9.em, flmaj, b), galax- 
ies in class b have the lowest Eddington ratio and radio 
power, and their nuclear and/or star formation activity 
even appears to be suppressed relative to the radio quiet 
( RQ) galaxies that have similar luminosities and mass 
(^ Table©. 

• Different researchers usually have adopted somewhat 
different definitions for the FR I/II types. As our pro- 
posed classification scheme is based on both radio mor- 
phology and nuclear emission line strength, and the 
original FR scheme is purely morphology-based, there 
is no one-to-one correspondence between the two. Nev- 
ertheless, given the similarities of class b and class a 
objects with r, < 0.8, and the large difference between 
these subsamples and the class a RGs with > 0.8 
OS.ll l, we can broadly identify the FR I type with the 
former a nd F R II with the latter sources (see the discus- 
sion in ^6.1l i. However, there is considerable overlap 
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in the distributions of physical properties for the three 
subsamples, and the trans ition f rom one FR type to the 
other is far from sharp (§ 3416.31 ). In particular, our find- 
ings do not support the previous claim that the two FR 
types occupy distinct regions in the radio luminosity- 
optical magnitude plane. 

• Although on average the ao.9.em objects are less massive 
than the other RG subsamples, they are still hosted by 
massive galaxies (MJ.' ' < M*). The distinction of this 
subset from the other RGs is far less dramatic than the 
blue-red bimodality of the general galaxy population. 

• To single out oo.g.em from the rest of the population 
in a statistically complete, low redshift RG sample for 
which optical emission line measurements are not avail- 
able, a possible approach is to select sources with 
ranked in the top 10% of the distribution. 

• Some of the objects in class b may form an evolutionary 
sequence, that is, they can be regarded as RGs seen at 
different stages of evolution, as evidenced by an anti- 
correlation between size and ( 36.21 Fig.[T3]l. A larger 
sample is needed to evaluate the effect of systematic 
uncertainties in the sample selection, however 

• Many different mechanisms must be at work for the 
generation of radio jets, but we suggest that the accre- 
tion rate onto the SMBH is the main driver for the dif- 
ferent radio morphologies, with host galaxy structure 
and/or density of the surrounding environment play- 
ing a secondary role. This is primarily motivated by 
the stark differences in the nuclear emission proper- 
ties of the RG subsets. At high accretion rates (e.g., 
[O III] A5007 line Eddington ratio > 10"^), the accretion 
mode is likely dominated by a geometrically thin, op- 
tically thick disk which could generate powerful, well- 
collimated jets that create strong hot spots at the edge 
of the lobes (i.e., the "classical double" morphology). 
At lower accretion rates, a radiatively inefficient accre- 
tion flow (RJAF) takes over; the outflows/jets from such 
accretion flows may not be as well-coUimated beyond 
the galactic nucleus scale, resulting in the "fat double" 
morphology (i.e., without obvious hot spots at the edge 
of the lobes). At slightly lower accretion rates (e.g., 
[O III] Eddington ratio < 10"^) and for galaxies with 
sufficiently high galactic density, a jet-dominated mor- 
phology is created ( 36.4b . 

• Based on the spectral properties of the galaxies, we sug- 
gest that outflows/jets from a RIAF may be more effi- 



cient in suppressing processes that cause star formation 
and/or nuclear activity than the jets from thin accretion 
disks. The latter could affe ct th e large scale surround- 
ings of the RGs, however ( 36.41 Fig.[8l). 

The advent of wide-field, uniform radio and optical surveys 
such as NVSS, FIRST, and SDSS makes it possible to pro- 
duce the large RG sample used here, and the classification 
scheme we propose. Although it is not clear if our scheme is 
more physically motivated than the existing ones (e.g., those 
of FR and OL89), our classification should be among the most 
objective and quantitative, and easily reproducible by other 
researchers. 

In this study we have only concerned ourselves with the 
extended sources with relatively "straight" lobes, that is, we 
have excluded wide-angle tail and narrow-angle tail objects. 
We have also left out the compact/point-like sources and radio 
quasars in the analysis. In a future publication we will com- 
pare the host properties of RGs of these other morphologies, 
which may provide further insights into the generation of the 
radio emission in galactic nuclei. 
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APPENDIX 

We noted in ^ that the measurements for both S and T may depend on the sensitivity and resolution of the radio data. We 
show below (Appendix lAli that the FIRST-based measurement of S should be robust for our RG sample. Our main results are 
also shown to be qualitatively insensitive to the data (FIRST or NVSS) from which total size T is measured (AppendixlBl). 

A. RESOLUTION OF FIRST IMAGES 

Part of the SDSS Stripe 82 has been observed with the VLA in the A configuration (1.5" FWHM) at 1.4 GHz, reaching rms 
0.07 mJy/beam. 15 RGs in our sample lie in the region covered by this deep survejQ (PI: G. Richards; Hodge et al. 2010, in 
preparation). Their redshifts range from 0.04 to 0.25, which is representative of our sample (Fig.|6l panel j). We compare the S 
and T measurements from these data with those from FIRST (hereafter with the subscripts V and F, respectively). Among the 15 



''^ Ihttp : / /www. physics ■ drexel ■ edu/ $\sim$gtr/vla/ stripe82/| 
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Fig. 15. — Median values of several physical properties as a function of rs (c.f. Fig.|6]and Fig.|2), with T derived from NVSS (magenta and cyan triangles) 
and those with T from FIRST (red and blue points). For class h RGs (blue and cyan points), we plot the results with negative values of rj to avoid cluttering the 
Figure. 



objects, five are in class b. For the ten class a RGs, the ratio Sy /Sf has a mean of 0.95 and a scatter of 0. 15. We also find that the 
total size measurements are very close, which suggests that the rj value derived from FIRST is robust against resolution issues. 

For four of the class a RGs, the HSB coincides with prominent jets in the host galaxy in the deep VLA images, and are 
classified as class b or c, depending on the central-to-total flux ratio (|3]). We note, however, that some of the flux from the lobes 
is resolved out in these A-array maps, and therefore ideally one should combine both the high and low resolution data to measure 
the proportion of flux that is in the jet component which only shows up in high resolution maps. These objects tend to have 
rj- values higher than the majority of the class a RGs, based on FIRST data; that is, FIRST does detect the central component. 
Because of the poorer angular resolution of FIRST, the jets are not as prominent as in the A-array maps. 

For the current analysis, we acknowledge the possibility that some of our objects which we have put into class a may in fact 
belong to class b if measured with better data. Using the distribution of (FlRST-based) r, and ry of the above four class a RGs, 
we estimate that 13% of class a objects may be subject to this misclassification. 

Although the redshift distribution for class a is quite consistent across the bins, class b RGs with higher r, are on average 
at slightly higher redshift (Fig. |6] panel j). Given the inherent sensitivity of the classification on the resolution, this is perhaps 
not surprising. For the five class b objects with deeper, higher-resolution VLA data, two have rj « 0.5 based on FIRST, and have 
redshift of 0.224 and 0.252, respectively. At three times finer resolution than FIRST, these RGs remain jet-dominated (i.e., class 
b), although their decreases. Given that the properties of class b objects do not vary much with respect to rj, we conclude that 
our results for class b should be robust (except for the possible addition of class a RGs with better measurements). 

B. MEASUREMENT OF TOTAL SIZE 

We have the option of using either NVSS or FIRST data to measure the total size T of the radio sources, and have chosen to use 
the latter as the default (except for the 41 cases where the FlRST-based sizes are much less than those from NVSS, presumably 
due to the insensitivity of FIRST to diffuse emission; see We have repeated our analysis with NVSS-based T measurements. 
If the fitted size of the major axis of radio sources from NVSS is only an upper limit, we exclude the sources from the sample; in 
addition, a minimal size of T = 50" (rather than T = 30" as adopted in is imposed, and therefore the sample size (797 RGs) 
with NVSS-based measurement is smaller. 

Variations of physical properties as a function of r,, analogous to Figures|6H8j are shown in Figures[T5]and[T6] for both FlRST- 
based and NVSS-based results. The main difference is in the total size Tp of the sources (Fig. [16] lower left panel): those derived 
from NVSS are larger than the ones based on FIRST, as expected. Since S is still measured using FIRST data and thus remains 
unchanged, the NVSS-based r, values are systematically lower than the FlRST-based ones; any trends with seen in Figure s|6]-[8] 
would therefore appear "stretched" horizontally and shifted towards low a bit. 
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TABLE 4 

Basic properties of Radio Galaxies 



SDSS ID 


VAGC ID" 


Class 


RA 


Dec 




log/'i.4 




cr 


''cff 






-f-OIII S0.5 








(J2000) 


(J2000) 




(W/Hz) 




(km/s) 


(kpc) 


(kpc) 






587731186743967908 


1659866 


b 


7.247278 


0.333494 


0.2222 


24.44 


-23.063 


331 


11.5 


111 


0.48 


■ ■ ■ 22.9 


588015510345613463 


1699168 


b 


5.794018 


0.976323 


0.2270 


24.68 


-23.039 


266 


9.2 


120 


0.24 


■■■ 15.6 


587724232636498080 


354077 


a 


7.205604 


14.979081 


0.0977 


24.96 


-21.627 


177 


6.8 


368 


0.56 


1.4 


587731185132503151 


1652016 


a 


5.281742 


-0.925413 


0.1082 


24.52 


-22.985 


268 


9.5 


117 


0.69 


14.2 



Note. — The full table is available and kept up-to-date at http://member.ipmu.jp/yen-ting.lin/RG/index.html 
NYU VAGC DR6 object ID. 



Using NVSS-based T measurements, we still find that class a objects with the highest r, stand out from the rest of the sample, 
although the division is now at rj w 0.7. The subtle difference in properties between class b and class a RGs with < 0.7 also 
remains. 

A possible concern of using either NVSS or FIRST to measure T is the relatively high SB limits of these surveys. One 
could imagine that a FD source with HSB spots far from the edge of the lobes would appear as high objects if observed with 
insufficient depth. We have checked against the NASA/IPAC Extragalactic Database (NED) to look for archival radio images for 
our ao.9 objects. Only 8 RGs 10%) have been imaged with decent data from the literature, and all of them would still have 
high Vs in those deeper maps. Even if such a bias due to the depth of the surveys exists, we suspect the difference in morphology 
between a bona fide hot spot at the edge of a lobe and a HSB region within a lobe for lower-r, RGs would be obvious enough in 
FIRST images s o tha t a visual inspection would be able to pick up such cases. 

We noted in ^6.11 that the combination of the presence of hot spots at the edge of the lobes and the high accretion rate as 
indicated by the presence of emission lines seems to be a pretty robust indicator for the type of RGs corresponding to FR 11, or the 
CDs. Since r, is unfortunately somewhat resolution-dependent, and emission line properties in the optical are not always easily 
available, perhaps a more objective approach to single out CDs from the rest of the population in a statistically complete, low 
redshift RG sample (so that our results are fully applicable) is to select sources with r, ranked in the top 10% of the distribution. 
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